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Abstract 



We present limits on the fraction of large angle microwave anisotropies which 
could come from tensor perturbations. We use the COBE results as well as 
smaller scale CMB observations, measurements of galaxy correlations, abun- 
dances of galaxy clusters, and Lyo absorption cloud statistics. Our aim is 
to provide conservative limits on the tensor-to-scalar ratio for standard in- 
flationary models. For power-law inflation, for example, we find T/S < 0.52 
at 95% confidence, with a similar constraint for <fP potentials. However, for 
models with tensor amplitude unrelated to the scalar spectral index it is still 
currently possible to have T/S > 1. 

PACS numbers: 98.80.-k, 98.70.Vc. 
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I. INTRODUCTION 



The presence of a primordial gravitational wave perturbation spectrum was an early 
prediction of inflationary models of the big bang [0 . However, it was not until the results of 
the COBE satellite mission that it became possible to begin to meaningfully constrain the 
tensor contribution to the overall perturbation spectrum 0-§f. In an early result, Salopek 
H found that, assuming power-law inflation, tensors must contribute less than about 50% 
of the cosmic microwave background (CMB) fluctuations at the 10° scale. 

Since that time, ground- and balloon-based experiments have begun to fill in the smaller- 
scale regions of the CMB power spectrum. These scales are crucial for constraining the 
gravity wave contribution because the tensor spectrum is expected to be negligible on scales 
finer than ~ 1°, and therefore large-scale power greater than that expected for scalars can be 
attributed to tensors. Markevich and Starobinsky have set some stringent limits on the 
tensor contribution. For example, they found that the ratio of tensor to scalar components 
of the CMB spectrum is T/S < 0.7 at 97.5% confidence for a flat, cosmological-constant-free 
universe with H = 50 kms _1 Mpc _1 . However, their analysis used a limited CMB data set 
and considered only a restricted set of values for the cosmological parameters. Recently 
Tegmark || has performed an analysis using a compilation of CMB data, and found a 68% 
upper confidence limit of 0.56 on the tensor to scalar ratio. In this work specific inflationary 
models were not considered, but a number of parameters were allowed to vary freely. In 
another recent study Lesgourgues et al. || analysed a particular broken-scale-invariance 
model of inflation with a steplike primordial perturbation spectrum, and found that the 
tensor to scalar ratio can reach unity. Melchiorri et al. [TIJ placed limits on tensors allowing 
for a blue scalar spectral index, and indeed found that blue spectra and a large tensor 
component are most consistent with CMB observations. 

Our aim here is to provide a more comprehensive answer (or set of answers) to the 
question: how big can T/S be? We present constraints on tensors for specific models of 
inflation as well as for freely varying parameters. In all cases we marginalize over the 
important, but as yet undetermined, cosmological parameters. We use both COBE and 
small-scale CMB data, as well as information about the matter power spectrum from galaxy 
correlation, cluster abundance, and Lyman a forest measurements. We refer to these various 
measurements of the power spectra as "data sets" . We additionally consider the effect, for 
each data set, of various observational constraints on the cosmological parameters, such 
as the age of the universe, cluster baryon density, and recent supernova measurements. 
We refer to these constraints as "parameter constraints" (this separation between "data 
sets" and "parameter constraints" is somewhat subjective, but dealt with consistently in 
our Bayesian approach; it is conceptually simpler to consider power spectrum constraints 
as measurements with some Gaussian error, while regarding allowed limits on cosmological 
parameters as restrictions on parameter space). Finally we consider what implications our 
results have for the direct detection of primordial gravity waves. 
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II. INFLATION MODELS 



Our goal is to provide limits on the tensor contribution to the primordial perturbation 
spectra using a variety of recent observations. In models of inflation, the scalar (density) 
and tensor (gravity wave) metric perturbations produced during inflation are specified by 
two spectral functions, A$(k) and A^(k), for wave number k. These spectra are deter- 



mined by the inflaton potential V(4>) and its derivatives [11]. However, when comparing 



model predictions with actual observations of the CMB, it is more useful to translate the 
inflationary spectra into the predicted multipole expansions of the CMB temperature field: 
AT/T(9, <p) = J2em a emYem{6,<j>), where Yi m (8,(f>) are the spherical harmonics. The spec- 
trum Ce = (\ae m \ 2 ) can be decomposed into scalar and tensor parts, Cg = Cf + Cj. In 
the literature the tensor to scalar ratio is conventionally specified either at I = 2 or in the 
spectral plateau at i ~ 10 — 20. Here we have chosen the i = 2 or quadrupole moments of 
the temperature field, and write S = 5Cf/47r and T = 5CJ/4it as usual ]T2"| . 

In order to constrain the tensor contribution T/S, we need to specify the particular 
model of inflation under consideration. This is because the model may provide a specific 
relationship between the ratio T/S and the scalar spectral index ns- Except in Sec. |VT| we 
only consider spatially flat inflation models (i.e. Qq + = 1, where Qq and Q\ = A/(3Hq) 
are the fractions of critical density due to matter and a cosmological constant, respectively). 
In addition, we do not consider the "quintessence" models |13|,14|, where a significant fraction 
of the critical density is currently in the form of a scalar field with equation-of-state different 
from that of matter, radiation, or cosmological constant (although it would not be difficult 
to extend our results for explicit models with recent epoch dynamical fields). We will also 
restrict ourselves to models which use the slow-roll approximation, and incorporate only a 



single dynamical field - a class of models sometimes called "chaotic inflation" [15]. This is 
not as restrictive as it might sound, since most viable inflationary models are of this form. 
Although some genuinely two-field models are known [16[], many multi-field models, the 



"hybrid" class, have only one field dynamically important and in these cases we effectively 
regain the single field case |T7|Jl8|| . In addition, theories which modify general relativity (e.g. 
"extended" inflation |I^,[2tJ) can often be recast as ordinary general relativity with a single 
effective scalar field pl"|,|i~8f . 

It is convenient to classify inflationary models as either "small-field" , "large-field" , or the 
already mentioned hybrid models |p2| . Small-field models are characterized by an inflaton 
field which rolls from a potential maximum towards a minimum at (cf)) ^ 0. These models 
generally produce negligible tensor contribution, but may result in the spectral index 



differing significantly from scale invariance [11]. In hybrid models, the important scalar field 



rolls towards a potential minimum with non-zero vacuum energy. These models also typically 
have very small T/S, and the scalar index can be greater than unity The large- field 
models involve so-called "chaotic" initial conditions, where an inflaton initially displaced 
from the potential minimum rolls towards the origin. Large-field models can produce large 
T/S and 1 — ns, and these are the models considered in this paper. This is not to say that 
small-field and hybrid models are not interesting; on the contrary, current views of inflation 



in the particle physics context suggest that T/S is expected to be small |23|| . However, 
large-field models must be considered when examining the observational evidence for a large 
tensor contribution. 
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It is also worth pointing out that we could construct models with a dip in the scalar 
power spectrum at large scales which compensates for the tensor contribution. Although we 
have not explored detailed models, we imagine that in principle models could be constructed 
with arbitrarily high T/S. We consider all such models with features at relevant scales to 
be unappealing unless there are separate physical arguments for them. 

In addition to considering models with free scalar index and tensor contribution T/S, we 
shall thus focus on two classes of inflationary models which can be considered representative 
of those predicting large gravity wave contributions. Both are restricted to "red" spectral 
tilts, ns < 1. The first, "power-law inflation" (PLI) pi||?5| , is characterized by exponential 
inflaton potentials of the form 



V(4>) oc exp 



l6lT<j) 2 



(1) 



and results in a scale factor growth a(t) oc t q , hence the name. For PLI the tensor-to-scalar 
ratio in fc-space can be calculated exactly as a function of ns: 



A\(k) _ 1 - n s 
IP) = 3 - n s 



(2) 



Note the tensor contribution is directly related to the scalar spectral index ns, which is 
further related to the tensor spectral index np = n$ — 1 in this model. Converting from 
fc-space to the observed anisotropy spectrum introduces a dependence on the cosmological 



constant which can be approximated by [12] 

T/S = -7n [0.97 + 0.58n + 0.25Q A - (l + l.ln + 0.28n 2 ) Q 



(3) 



where n = ns — 1 = The dependence on A arises because of different evolution for 
scalars and tensors when A dominates at late times. The dependence on other cosmological 



parameters is negligible [|T^ . 

We also consider the large-field polynomial potentials, 



for integral p > 1 



V{4>) oc <jf, (4) 

. In this case both ns and np are determined by the exponent p JTT 

2p + A 



1 - 



p + 200' 



(5) 



rip 



-2p 
p + 200' 



(6) 



The tensor index may be related to T/S through the consistency relation [12 



T 

S 



A 0) 

7 j ± 

is 



(7) 
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where the cosmo logical parameter dependence, again dominated by , can be approximated 
by 

ffi = 1.04 - 0.82ft A + 2fi A , (8) 

/4 0) = 1.0 - o.o3^ A - o.m 2 A . (9) 

As a third possibility, we will also consider models with scalar index varying over the 
range ns = 0.8 — 1.2, but with an independently varying tensor contribution T/S. 



III. MICROWAVE BACKGROUND ANISOTROPIES 



In order to evaluate likelihoods and confidence limits for T/S based on CMB measure- 
ments, we performed \ 2 fits °f model Ci spectra to CMB data. We did this for a set of 
"band-power" estimates of anisotropy at different scales, and separately for the COBE data 
themselves. For our first approach, we used a collection of binned data to represent the 
anisotropies as a function of t. Specifically we took the flat-spectrum effective quadrupole 
values listed in Smoot and Scott |27j and binned them into nine intervals separated logarith- 



mically in t. We chose this simplified approach since we anticipated a large computational 
effort in covering a reasonably large parameter space. The use of binned data has been 
shown elsewhere [28] to give similar results to more thorough methods. If anything, there is 



a bias towards lowering the height of any acoustic peak, inherent in the simplifying assump- 
tion of symmetric Gaussian error bars |[29|| ; for placing upper limits on T/S our approach is 



therefore conservative. We are also erring on the side of caution by using the binned data 
only up to the first acoustic peak, neglecting constraints from detections and upper limits 
at smaller angular scales. 

We ignored the effect of reionization on the Cg spectra. Reionization to optical depth r 
reduces the power of small-scale anisotropies by e~ 2r . Thus, in placing upper limits on T/S, 
it is conservative to set r = 0. 

A fitting function for the spectrum, valid up to the first peak at I ~ 220, has been 



provided by White 30 




c i(y) = \7z) c l {v = o), (io) 



where v is the (nearly) degenerate combination of cosmological parameters 

v = n s - 1 - 0.32 ln(l + 0.76r) + 6.8(Q B h 2 - 0.0125) - 0.371n(2/i) - 0.161n(fi )- (H) 

Here r = lACj Q /Cf Q is the tensor to scalar ratio at i = 10, normalized to provide r = T/S 
for Q = 1 and ns — > 1. The parameter h is defined through H = 100/i kms _1 Mpc _1 , 
and Ob is the fraction of the critical density in baryons. Thus the standard CDM (sCDM) 
spectrum is specified by v = 0. We found that the fl A dependence of r can be well captured 
by introducing the rescaled variable r', defined by 

r ' = 0.94 + 1.105fi 3 A ' 75 ' (12) 



5 



and setting r' = T / S. 

We fitted the model spectra of Eq. fllTf ) to the binned data as follows. For each com- 
bination of parameters (h, ^h 2 , £Iq, ns, T/S) we normalized the model spectrum to the 
binned data, and evaluated the likelihood C(h,fl B h 2 ,fl ,ns,T/S) oc exp(— x 2 /2). Next this 
likelihood was integrated, uniformly in the parameter, over the ranges of h = 0.5 — 0.8, 
n B h 2 = 0.007 - 0.024, and fi = 0.25 - 1, subject to the constraints of Eq. (§) for PLI and 
Eqs. (|5|), (^), and (|7|) for polynomial potentials. For the case of free T/S, the scalar index 
was varied in the range n$ = 0.8 — 1.2. Finally the resultant C(T/S) was normalized to a 
peak value of unity and the 95% confidence limits evaluated. We tried to choose reason- 
able ranges for the prior probability distributions of the "nuisance parameters" , guided by 
the current weight of evidence. We checked that mild departures from our adopted ranges 
lead to only small modifications to our results. However, we caution that our conclusions 
will not necessarily be applicable for models which lie significantly outside the parameter 
space we considered. In addition, note that according to Eq. flTT|) we can crudely estimate 
an upper limit on T/S by combining the observational lower limit on v with the maximal 
baryon density and minimal and h from our parameter ranges. However, this turns out 
to be an overly conservative estimate: for example, for ns = 1, and using a lower limit of 
v = —0.2, Eq. flTTp gives an upper limit of T/S = 3.8, compared with the limit T/S = 1.6 
from Sec. |V11| . 

For the separate constraint from the COBE data, we used the software package CMB- 
FAST jH] to calculate likelihoods based only on the COBE results at large scales. CMB- 
FAST calculates the spectrum using a line-of-sight integration technique. It then calculates 
likelihoods by finding a quadratic approximation to the large scale spectrum and using the 
COBE fits of Bunn and White [j32fl . These likelihoods were integrated and 95% limits calcu- 
lated as above, except that the baryon density was fixed at Q B = 0.05 to save computation 
time (and since Qb has negligible affect at these scales). The results of this procedure are 
presented in Sec. |V11. 



IV. LARGE-SCALE STRUCTURE 



A. Galaxy correlations 



We next applied observations of galaxy correlations to constrain T/S indirectly through 
the power spectrum of the density fluctuations, A 2 (A;). The power spectrum A 2 (k) is ex- 
pressed, following Bunn and White |[32| , by 

f rh\ 3+ns 

A 2 W=4Kr T2 (k). (13) 



Here 5n is the (O , ns, and r dependent) normalization described in Sec. [IV B| , and T(k) is 
the transfer function which describes the evolution of the spectrum from its primordial form 
to the present. 

We explicitly used for the transfer function the fit of Bardeen et al. |33| . 
ln(l + 2.34g) N , _ , a , 



T(q) 



2.Uq 



1 + 3.89g + (16. IqY + (5.46gf + (6.71g) 4 



(14) 
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with the scaling of Sugiyama 



P4 



MT 70 /2.7K) 2 



fio/i 2 exp 

(-n B - ^/V - 5 ^b/^o) ' 

Here T 70 is the temperature of the CMB radiation today. 

We performed \ 2 fits of the (unnormalized) model power spectrum given by Eqs. ([13]) - 
(f5|) to the compilation of data provided in Table I of Peacock and Dodds []3~5| |, excluding 



their four smallest scale data points. These points were omitted because, while there are 



theoretical reasons |36[ to expect that the galaxy bias approaches a constant on large scales, 
at the smallest scales the assumption of a linear bias appears to break down |J7|. Here we 
are fitting for the shape of the matter power spectrum, ignoring the overall amplitude, since 
the normalization is complicated by the ambiguities of galaxy biasing. 



The fitting was performed in exactly the same way as was described in Sec. |T| for the 
binned microwave anisotropics. Namely the model curves were normalized to the Peacock 
and Dodds data, the integrated likelihood was calculated, and the 95% confidence limits for 
n s were evaluated. Since the shape of the power spectrum [Eq. (jl5|)] is independent of the 
tensor amplitude, this technique can only provide limits on T/S when T/S is determined 
by the spectral index. That is, the galaxy correlation data can only constrain T/S for our 
PLI and cjf cases, using the relationships [Eq. (|]) or Eqs. (Q), (||), and (|7|)] between T/S and 
n s . 



B. Cluster abundance 



A very useful quantity for constraining the amplitude of the power spectrum is the 
dispersion of the density field smoothed on a scale R, defined by 



a 2 {R) 



W 2 {kR)A 2 {k) 



dk 



(16) 



Here W{kR) is the smoothing function, which we take to be a spherical top-hat specified by 

sin(kR) cos(kR) 



W(kR) 



(17) 



(kRf (kR) 2 

Traditionally the dispersion is quoted at the scale 8h~ l Mpc, and given the symbol a 8 . 



For our experimental value we used the result of Viana and Liddle |38| , who analysed the 
abundance of large galaxy clusters to obtain 



0.56Qq 0A7 , 



(18) 



with relative 95% confidence limits of — 18 f^' 21ogl0 ^° and +20 21ogl0 percent. Several 
other estimates have been published; the one we used is fairly representative, and with a 
more conservative error bar than most. 

To compare this experimental result with the model value predicted by Eq. (|i~6"D, we must 
fix the normalization <5 H - We used the result of Liddle et al. who fitted <5h using the 
COBE large scale normalization to obtain 
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10 5 S H (n s ,n ) = 1.94fi -°- 785 - 051n ^exp[/(n s )], (19) 

where 

r/ \ [ — 0.95n — 0.169n 2 , No tensors, , nnS 

/(ns) = \l.00n + 1.97n 2 , PLI. (20) 

For the case of non-PLI tensors, we used the fitting form of Bunn et al. 

10 5 5n = 1.91 o-o-so-o-osmno exp(-l.Oln) + Q _ Q Q3 

^1 + (0.75 - 0.13fii) r 

We calculated likelihoods for our model ag using a Gaussian with peak and 95% limits 
specified by Eq. (|T8|), and then integrated £ and found limits for T/S as in the binned 
microwave case. 



C. Lyman a absorption cloud statistics 

Another measure of the amplitude of the matter power spectrum has been obtained 



recently by Croft et al. [fil l, who analysed the Lyman a (Lya) absorption forest in the 
spectra of quasars at redshifts z ~ 2.5. These results apply at smaller comoving scales than 
the cluster abundance ag measurements, and hence are potentially more constraining. Croft 
et al. found 

A 2 (A; P ) = 0.57^ (22) 

at la confidence, where the effective wavenumber k p = 0.008(km s -1 ) -1 at z = 2.5. 

These results cannot be directly compared with the model predictions of Eq. fljjD, because 
Eq. fll3|) provides its predictions for the current time, i.e. z = 0. To translate to z = 2.5, 
we must first convert the model k from the comoving Mpc -1 units conventionally used in 
discussions of the matter power spectrum to (km s -1 ) -1 at z = 2.5, using 

^(kms-^H^MMpc- 1 ], (23) 

where 



H(z) =# \M)(1 + ^) 3 + ^a (24) 

for flat universes. 

Next, we must consider the growth of the perturbations themselves. In a critical density 
universe (and assuming linear theory), the growth law is simply A 2 (k, z) = A 2 (k, 0)(l + z)~ 2 . 
As fl\ increases, the growth is suppressed, and this can be accounted for by writing 

A'(*.») = A' (fc ,0) ™ 1 ■ (25) 
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where the growth suppression factor g(fl) can be accurately parametrized by ||2 



_ 5/1 209ft Q 2 

o(ft) = -ft — + + ft 4/7 , (26) 

yv ; 2 \70 140 140 y ' V ; 

and the redshift dependence of ft is given by 

n M = *T=n£$W (27) 

all for spatially flat universes. 

We calculated likelihoods using the normalized model predictions of Eq. flHf), translated 

to z = 2.5 as described above, and then obtained limits for T/S as in the cluster abundance 
case. 



V. PARAMETER CONSTRAINTS 
A. Age of the universe 



In flat A models, the age of the universe is j|2 



tu _ 2 sinh- 1 (y/^My (2s) 
3H v / ft~A 

During the integration of the likelihoods, we investigated the effect of imposing a constraint 
on the parameters h and ft , so that regions of parameter space corresponding to ages below 
various limits were excluded. This simply corresponds to a more complex form for the priors 
on the parameters. The precise limit on the age of the universe is a matter of on-going debate 



(e.g. [f43| , [44| ) . A lower limit of around 11 Gyr now seems to be the norm, so we considered 
this case explicitly. We also considered the effect of a more constraining limit of 13 Gyr, 
still preferred by some authors. 



B. Baryons in clusters 

Recent measurements of the baryon density in clusters have suggested low fto for con- 
sistency with nucleosythesis. We chose to use the results of White and Fabian [|^j for the 
baryon density 

^ = (0.056 ± 0.014)/T 3/2 , (29) 
fto 

where the errors are at the la level. We explored the implications of applying this constraint 
during the likelihood integrations, by adding a term 

1 2 (30) 



/ /i 3 ^ B M- 0-056 ^ 
V 0.014 t 



to each value of % 2 . 
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C. Supernova constraints 



Measurements of high-z Type-la supernovae (SNe la) are in principle well-suited to con- 
straining Q on the assumption of a flat A universe, since such measurements are sensitive to 
(roughly) the difference between Q and We used the experimental results of Filippenko 
and Riess of the High-z Supernova Search team p6], who found for flat A models 



Q = 0.25 ±0.15 (31) 
at la confidence. We also investigated the effect of applying this constraint as above. 



VI. OPEN MODELS 

For models with open geometry the situation is more complicated, and so we restrict 
ourselves to a brief discussion here. In addition to the added technical complexity involved 
in working in hyperbolic spaces, the presence of an additional scale, the curvature scale, 
renders ambiguous the meaning of scale-invariant fluctuations. For the most obvious scale- 
invariant spectrum of gravity wave modes, the quadrupole anisotropy actually diverges! For 
this reason one requires a definite calculation of the fluctuation spectrum from a well realized 



open model. The advent of open inflationary models [[47] has allowed, for the first time, a 
calculation of the spectrum of primordial fluctuations in an open universe. As with all 
inflationary models, a nearly scale-invariant spectrum of gravitational waves (tensor modes) 



is produced [ 48|J49 |. The size of these modes in fc-space, and their relation to the spectral 
index, is not dissimilar to the flat space models we have been considering. In the inflationary 
open universe models the spectrum of perturbations is cut-off at large spatial scales, leading 
to a finite gravity wave spectrum. However, the exact scale of the cutoff depends on details 
of the model, introducing further model dependence into the £-space predictions. 

Since gravitational waves provide anisotropies but no density fluctuations, their presence 
will in general lower the normalization of the matter power spectrum (for a fixed large angle 
CMB normalization). Open models already have quite a low normalization |5Tl,5"i|, so the 
most conservative limits on gravity waves come from models which produce the minimal 
tensor anisotropies, i.e. where the cutoff operates as efficiently as possible. The COBE 



normalization for such models with PLI is [5"2~| 



10 5 5 n = 1.95 fi Q 0-35-0.19 In n +0.15n exp A Q2 ~ + UQjjfj . ( 32 ) 

Combining this normalization with the cluster abundance gives a strong constraint on T/S. 
We show in Fig. |l]the 95% CL upper limit on T/S as a function of Qq i n these models. 



VII. RESULTS 

Figure |2] presents likelihoods, integrated over the parameter ranges described above, and 
plotted versus T/S, for the various data sets, and specifically for PLI models. For the curve 
labelled "combined" , likelihoods for each data set (except the COBE data) were multiplied 
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FIG. 1. 95% upper confidence limits on T/S for the open models described in the text. The 
cluster abundance data set was used with no parameter constraints. 

together before integration. (Including the COBE data would have been redundant, since 
the binned CMB set already contains the COBE results.) Thus the "combined" values 
represent joint likelihoods for the relevant data sets, on the assumption of independent 
data. Note that the combined data curve of Fig. differs significantly from the product of 
the already marginalized curves for the different data sets, which indicates that parameter 
covariance is important here. Also, the maximum joint likelihood in Fig. |2| corresponds to 
X 2 — 9, which indicates a good fit for the 15 degrees of freedom involved. 

Figure |3| displays integrated likelihoods versus T / S for each data set and for the combined 
data, again on the assumption of PLI. The effect of each parameter constraint is illustrated. 
The COBE data shape constraint is very weak, and exhibits essentially no cosmological 
parameter dependence, as expected. Thus the parameter constraints have little effect on 
the likelihoods, and the curves are not shown here. Cluster abundance is not much more 
constraining than the COBE shape, but exhibits considerably stronger cosmological param- 
eter dependence, and hence is affected substantially by the various parameter constraints. 
The matter power spectrum shape constraint is so weak that we do not plot it here. The 
strongest constraint comes from the binned CMB data, and indeed these data dominate the 
joint results. 

We can understand the general features of the large parameter dependence exhibitted by 
the likelihoods for the matter spectrum data sets as follows. Near sCDM parameter values, 
it is well known that the matter power spectrum contains too much small-scale power when 
COBE-normalized at large scales. The presence of tensors improves the fit at small scales 
by decreasing the scalar normalization at COBE scales. Reducing /i or O , however, also 
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1 1.5 

T/S 

FIG. 2. Integrated likelihoods versus T/S for the various data sets: short-dashed, long-dashed, 
dotted, dash-dotted, and solid curves represent respectively COBE, binned CMB, cluster abun- 
dance, Lya absorption, and combined data. All curves are for PLI inflation, using the priors 
discussed in Sec. [II, with no additional parameter constraints. 
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FIG. 3. Integrated likelihoods versus T/S for PLI and for the various data sets; clockwise 
from upper left: binned CMB, cluster abundance, Lya absorption, and combined data. Solid, 
short-dashed, long-dashed, dotted, and dot-dashed curves represent no constraint, to > H Gyr, 
to > 13 Gyr, cluster baryon fraction, and SNe la parameter constraints, respectively. 
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TABLE I. 95% confidence limits on T/S for various data sets and parameter constraints, and 
for power-law inflation. "2.5+" represents no constraint (for values as high as this we do not expect 



our approximations to be adequate in any case). 

Data set No constraint to > H Gyr to > 13 Gyr Cluster baryon Supernova 

COBE 21 21 21 21 21 

binned CMB 0.60 0.62 0.67 0.63 0.65 

galaxy correlations 2.5+ 2.5+ 2.5+ 2.5+ 2.5+ 

cluster abundance 1.8 1.4 1.1 1.5 1.1 

Ly a absorption 1.6 1.1 0.82 1.3 0.96 

combined 0.52 0.52 0.51 0.53 0.47 



TABLE II. 95% confidence limits for T/S as for Table | but for <jjP inflation. 



Data set No constraint to > H Gyr to > 13 Gyr Cluster baryon Supernova 
binned CMB 063 064 067 064 065 
galaxy correlations 2.5+ 2.5+ 2.5+ 2.5+ 2.5+ 
cluster abundance 2.1 1.5 1.2 1.8 1.2 
Lya absorption 1.8 1.2 0.87 1.5 1.1 
combined 0.49 0.49 0.49 0.50 0.45 



decreases the power at small scales, improving the fit over sCDM, and thus reducing the 
need for tensors. When an age constraint is applied, we force the model towards lower h 
and Q according to Eq. (p8|), and hence towards lower T/S, as is seen in Fig. |[ The cluster 
baryon and supernova constraints similarly move us to smaller Q . 

Figure ^ displays likelihoods versus p for <ff inflation, while Fig. |5] presents likelihoods 
versus T/S for the case of free tensor contribution and ng = 1. In all plots, curves have been 
omitted for the very weakly constraining data sets. The curves of Fig. [| closely resemble 
those of Fig. [3|. This is because, for p>2, Eqs. and (0) give T/S ~ — 6.85n, which 

is similar to the PLI result of Eq. (|J). 

In Fig. H we see that the data are considerably less constraining, compared with the 
PLI case, when we allow T/S to vary freely. This was expected, since in the PLI case, 
the lowering of ng tends to enhance the effect of increasing T/S. We also expect that a 
blue scalar tilt would oppose the effect of tensors on the spectrum and hence allow larger 
T/S. Figure || illustrates this effect by plotting the 95% upper confidence limits on T/S 
versus scalar index with T/S free. Note that we cannot meaningfully constrain T/S here 
by marginalizing over ng, since the best fits to the spectrum remain good even for very blue 
tilts and very large T/S. 

Our confidence limits are summarized in Table | for the case of power-law inflation, 
Table || for polymomial potentials, and Table |T| for free T/S and ng = 1. In all cases 
95% upper limits on T/S are presented, after integrating over the ranges of parameter space 
specified in Sec. |II[ The row gives the data set used, while the column specifies the type of 
parameter constraint applied, if any. 
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10 20 30 20 40 60 

P P 

FIG. 4. Integrated likelihoods versus p for cfrP inflation and for various data sets; clockwise 
from upper left: binned CMB, cluster abundance, Lya absorption, and combined data. Solid, 
short-dashed, long-dashed, dotted, and dot-dashed curves represent no constraint, to > H Gyr, 
to > 13 Gyr, cluster baryon fraction, and SNe la parameter constraints, respectively. 
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FIG. 5. Integrated likelihoods versus T/S for T/S free and n$ = 1, for binned CMB (left) and 
combined data (right). Solid, short-dashed, long-dashed, dotted, and dot-dashed curves represent 
no constraint, to > 11 Gyr, to > 13 Gyr, cluster baryon fraction, and SNe la parameter constraints, 
respectively. 
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FIG. 6. 95% confidence limits on T/S versus ns, for T/S free and for binned CMB (dashed) 
and combined data (solid). No parameter constraints have been applied. 



TABLE III. 95% confidence limits on T/S as for Table | but for ns = 1 and T/S free. 



Data set No constraint to > 11 Gyr to > 13 Gyr Cluster baryon Supernova 
COBE 2^5+ 2^5+ 2^5+ 2^5+ 2.5+ 
binned CMB 1.6 1.6 1.8 1.6 1.8 
cluster abundance 2.5+ 2.5+ 2.5+ 2.5+ 2.5+ 
Ly a absorption 2.5+ 2.5+ 2.5+ 2.5+ 2.5+ 
combined 1.3 1.3 1.3 1.3 1.0 
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VIII. THE FUTURE 



The discovery of a nearly scale-invariant spectrum of long wavelength gravity waves 
would be tremendously illuminating. Inflation is the only known mechanism for producing 
an almost scale-invariant spectrum of adiabatic scalar fluctuations, a prediction which is 
slowly gaining observational support. In the simplest, "toy" , models of inflation a potentially 
large amplitude almost scale-invariant spectrum of gravity waves is also predicted. For 
monomial inflation models within the slow-roll approximation, detailed characterization of 
this spectrum could in principle allow a reconstruction of the inflaton potential [[11]]. This 



surely is a window onto physics at higher energies than have ever been probed before. 

Inflation models based on particle physics, rather than "toy" potentials, predict a very 



small tensor spectrum [E3|. However, essentially nothing is known about particle physics 
above the electroweak scale, and extrapolations of our current ideas to arbitrarily high 
energies could easily miss the mark. We must be guided then by observations. We have 
argued that observational support for a large gravity wave component is weak. Indeed 
observations definitely require the tensor anisotropy to be subdominant for large angle CMB 
anisotropies. One the other hand, it is still possible to have T/S ~ 0.5, and since it would 
be so exciting to discover any tensor signal at all we are led to ask: how small can a tensor 
component be and still be detectable? What are the best ways to look for a tensor signal? 



A. Direct detection 



The feasibility of the direct detection of inflation-produced gravitational waves has been 
addressed by a number of authors P,|53|-p^|, with pessimism expressed by most. 

The ground-based laser interferometers LIGO and VIRGO J59| will operate in the 
/ ~ 100 Hz frequency band, while the European Space Agency's proposed space-based 
interferometer LISA |6(| would operate in the / ~ 10~ 4 Hz band. Millisecond pulsar tim- 
ing is sensitive to waves with periods on the order of the observation time, i.e. frequencies 
/ ~ 10~ 7 — 10~ 9 Hz [j59 |. These instruments probe regions of the tensor perturbation spec- 
trum which entered during the radiation dominated era. Expressions for the fraction of 
the critical density due to gravity waves per logarithmic frequency interval can be found in 
|5^-p7|. Assuming that f2o = 1 in a PLI model, with the only relativistic particles being 
photons and 3 neutrino species, and taking the COBE quadrupole Q = T + S ~ 4.4 x 10~ n , 
one finds \E% 



GW (/)^ 2 = 5.1 x 10 



-15 



n T - 1/7 



exp 



+ -N 2 (dn T /d Ink) 



(33) 



where N = ln(k/Ho) and = —{T/S) /I is the tensor spectral index. 

Using Eq. (|33|) Turner |57| found that the local energy density in gravity waves is max- 
imized at T/S — 0.18 for / ~ 10~ 4 Hz. At this maximum, the local energy density is in 
the range Qgw^ 2 — 10~ 15 — 10~ 16 , which lies a couple of orders of magnitude below the 
expected sensitivity of LISA, and several orders below that of LIGO/VIRGO This is 
also well below the current upper limit of f2cw^ 2 < 6 x 10~ 8 (at 95% confidence) from pulsar 
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timing fl6l] . As T/S increases above 0.18, Qgw(/ ~ 10~ 4 Hz)/i 2 begins to decrease due to 
the increasing magnitude of the tensor spectral index. 

Recall that our joint data constraint for PLI gives T/S < 0.5, so our results predict that 
the inflationary spectrum of gravity waves from PLI is not amenable to direct detection. 



B. Limits from the CMB 

With the advent of MAP and especially the Planck Surveyor, with its higher sensitivity, 
detailed maps of the CMB are just around the corner. What do we expect will be possible 
from these missions? This question has been dealt with extensively before. Assuming a 
cosmic variance limited experiment capable of determining only the anisotropy in the CMB 
but with all other parameters known, one can measure T/S only if it is larger than about 



10% [62]. A more realistic assessment for MAP and Planck suggests this limit is rarely 
reached in practice [E3,E3]. 

However the ability to measure linear polarization in the CMB anisotropy offers the 
prospect of improving the sensitivity to tensor modes (for a recent review of polarization see 
|65fl). In addition to the temperature anisotropy, two components of the linear polarization 



can be measured. It is convenient to split the polarization into parity even (E'-mode) and 
parity odd (B-mode) combinations - named after the familiar parity transformation prop- 
erties of the electric and magnetic fields, but not to be confused with the E and B fields of 
the electromagnetic radiation. 

Polarization offers two advantages over the temperature. First, with more observables the 
error bars on parameters are tightened. In addition the polarization breaks the degeneracy 
between reionization and a tensor component, allowing extraction of smaller levels of signal 
Model dependent constraints on a tensor mode as low as 1% appear to be possible 



with the Planck satellite ||6^|j6^ , |67|j68[| . Extensive observations of patches of the sky from 
the ground (or satellites even further into the future) could in principle push the sensitivity 
even deeper. 

There is a further handle on the tensor signal however. Since scalar modes have no 
"handedness" they generate only parity even, or E-mode polarization p6| , |69| . A definitive 



detection of 5-mode polarization would thus indicate the presence of other modes, with 
tensors being more likely since vector modes decay cosmologically. Moreover a comparison 
of the -B-mode, .E-mode and temperature signals can definitively distinguish tensors from 
other sources of perturbation (e.g. ]7CT|). 

Unfortunately the detection of a I?-mode polarization will prove a formidable experi- 
mental challenge. The level of the signal, shown in Fig. ^ for T/S = 0.01, 0.1 and 1.0, 
is very small. As an indicative number, with T/S = 0.5, our upper limit, the total rms 
5-mode signal, integrated over £, is 0.24 /iK in a critical density universe. These sensitivity 
requirements, coupled with our current poor state of knowledge of the relevant polarized 
foregrounds make it seem unlikely a S-mode signal will be detected in the near future. 
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FIG. 7. i?-mode tensor polarization signal Cf^e. f° r T/S = 0.01, 0.1, and 1.0, with the remaining 
parameters specified as standard CDM. 



We have examined the current experimental limits on the tensor-to-scalar ratio. Using 
the COBE results, as well as small-scale CMB observations, and measurements of galaxy 
correlations, cluster abundances, and Lya absorption we have obtained conservative limits 
on the tensor fraction for some specific inflationary models. Importantly, we have considered 
models with a wide range of cosmological parameters, rather than fixing the values of fl , 
H , etc. For power-law inflation, for example, we find that T/S < 0.52 at the 95% confidence 
level. Similar constraints apply to (ff inflaton models, corresponding to approximately p < 
8. Much of this constraint on the tensor-to-scalar ratio comes from the relation between 
T/S and the scalar spectral index n s in these theories. For models with tensor amplitude 
unrelated to scalar spectral index it is still possible to have T/S > 1. Currently the tightest 
constraint is provided by the combined CMB data sets. Since the quality of such data are 
expected to improve dramatically in the near future, we expect much tighter constraints (or 
more interestingly a real detection) in the coming years. 
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